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PREFACE 


The University of Texas at San Antonio (UTSA) is a relatively new university. It 
was established in 1969 and opened for classes in 1973. As the only comprehensive public 
university serving the nation's ninth largest city, it was and is vital to San Antonio and the 
entire South Texas Region. In 1983, just seven years ago, an undergraduate engineering 
program was established at UTSA with the support of the community and its leaders. 
Today, all three undergraduate engineering programs are ABET accredited and serve about 
800 students, a significant percentage of whom are Hispanic. The future includes a new 
engineering building, providing new laboratory facilities and equipment, together with 
offices and laboratories, planned to open in January, 1991. Furthermore, a graduate 
program is planned at both M.S. and Ph.D. levels, and it is hoped that the first Master’s 
Degree students will be able to enroll in Fall, 1989. 

Naturally, the engineering research environment is just developing at UTSA. Now, 
thanks in great measure to the UT System support and this ongoing NASA grant, good 
progress is being made. Specifically, the purchase of a UT System CRAY-XM/P in 
March, 1986 and a second one in December, 1988 has provided a world-class analytical 
and numerical research environment not ordinarily available to a new university. As a 
result the UTSA Supercomputer Network Research Facility (SNRF) was developed by the 
principal investigator, Dr. Lola Boyce. This has allowed the successful completion of this 
research project, the first of its kind at UTSA. 

This NASA research grant has allowed two Mechanical Engineering students, 
Thomas Lovelace and Callie Scheidt, to work directly with the principal investigator. Dr. 
Boyce, providing them with a quality research experience they would otherwise probably 
not have had. Both students have expressed an interest in continuing their educations at the 
graduate level. 

In conclusion, and in view of the significant accomplishments in fundamental 
research, enhancement of the engineering research environment at UTSA, and direct 
support of Mechanical Engineering students, it is hoped that the proposed extension of this 
grant will receive favorable consideration at NASA. The principal investigator sincerely 
thanks NASA for funding this first year grant 


ABSTRACT 


This report presents the results of the first year of effort of a program of research 
conducted for NASA-LeRC by The University of Texas at San Antonio (UTS A). The 
research included development of methodology that provides a probabilistic treatment of 
lifetime prediction of structural components of aerospace propulsion systems subjected to 
fatigue. Material strength degradation models, based on primitive variables, include both a 
fatigue strength reduction model and a fatigue crack growth model. Linear elastic fracture 
mechanics is utilized in the latter model. Probabilistic analysis is based on simulation, and 
both maximum entropy and maximum penalized likelihood methods are used for the 
generation of probability density functions. The resulting constitutive relationships are 
included in several computer programs, RANDOM2, RANDOM3 AND RANDOM4. 
These programs determine the random lifetime, of an engine component, in mechanical 
load cycles, to reach a critical fatigue strength or crack size. The material considered was a 
cast nickel base-superalloy, one typical of those used in the Space Shuttle Main Engine 
(SSME). 
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1.0 INTRODUCTION 


This report presents the results of the first year effort of a research program entitled 
"Development of Advanced Methodologies for Probabilistic Constitutive Relationships of 
Material Strength Models." This research is sponsored by the National Aeronautics and 
Space Administration-Lewis Research Center (NASA-LeRC). The principal investigator is 
Dr. Lola Boyce, Associate Professor of Mechanical Engineering, The University of Texas 
at San Antonio (UTSA). The objective of the research program is the development of 
methodology that provides a probabilistic treatment of lifetime prediction of structural 
components of aerospace propulsion systems subjected to fatigue. 

Two material strength degradation models, based on primitive variables were 
developed as part of this first year effort: a fatigue crack growth model and a fatigue 
strength reduction model. The former model utilizes principles of linear elastic fracture 
mechanics while the latter is, recently developed at NASA-LeRC, quantifies the reduction 
of strength under cyclic loading, including elevated temperature treatment Probabilistic 
analysis is based on simulation, and both maximum entropy and maximum penalized 
likelihood methods are used for the generation of probability density functions that predict 
the random lifetime of a material typical of those used in the Space Shuttle Main Engine 
(SSME), namely a cast nickel base-superalloy. 

The resulting constitutive relationships are included in several computer programs, 
RANDOM2, RANDOM3, and RANDOM4. The programs were developed using both the 
NASA-LeRC and UTSA Supercomputer Network Research Facility (SNRF) Cray X-MP. 
New versions of the program accompany this report (see enclosed floppy disk), utilizing 
the new IMSL Ver. 10 subroutines. Thus, these new versions of the programs will 
execute on the current NASA-LeRC supercomputer facilities. Also the floppy disk 
contains sample problems to verify program performance at NASA-LeRC. 

Finally, a sensitivity study was carried out for the fatigue strength reduction model 
for the case of a relatively high mean stress and a relatively low constant amplitude 
alternating stress at failure. In addition to varying the stresses, the effect of temperature 
was also considered. A paper was produced documenting much of the effort of this first 
year research program. This paper is entitled "Probabilistic Constitutive Relationships for 
Cyclic Material Strength Models", by L. Boyce and C.C. Chamis. It was presented at the 
29th Structures, Structural Dynamics and Materials Conference, Williamsburg, VA, April, 
1988 and is published in the Proceedings. It has also been submitted to the AIAA Journal 
of Propulsion and Power. 
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2.0 FATIGUE CRACK GROWTH MODEL 


2.1 Background 

Fatigue crack growth data are usually presented as cycles, N, to reach a particular 
crack length, a. The initial crack size is ai. It is generally accepted that under constant 
amplitude alternating stress, fatigue crack growth can be related to stress intensity through a 
first order differential equation 1 

da/dN = C(AK) m (1) 

where C is a material parameter, m is a material property (often a constant) and AK is the 
stress intensity range. Stress intensity range is given by 

AK = YAo/ia 

where Y is a constant dependent upon component and crack geometry and Act is the 
constant amplitude alternating stress. Therefore, equation (1) can be written as 

da/dN = C(YACT/ia) m 


or. 


da/dN = C Y m Ao m 7t m/2 a m/2 . (2) 

Equation (2) can be integrated, from the initial crack length, ai, to the final crack 
length, af, to yield N, the number of cycles. The result is 


CY m 7t m/2 Ao m 

Thus, equation (3) gives the "cycles to reach a given crack length." 

Metallurgical evidence indicates that casting pores play a significant role in the high- 
cycle fatigue life of cast nickel base-superalloys, especially at high temperatures. 2 The 
location and size of these fatigue crack-initiating pores vary greatly from one aerospace 
propulsion system component to another. This accounts for the large variability in fatigue 
life and leads to consideration of fatigue crack growth as a random phenomenon. 

Fatigue life directly relates to casting pore size, and pore size can be used to 
determine initial crack size, aj. Thus, utilizing principles of both probabilistic analysis and 
fatigue crack growth, a quantative probabilistic constitutive relationship between fatigue life 
and fracture mechanics parameters can be developed. Using the "randomized equation" 
approach, the fatigue crack growth model, given by equation (3) has the following form: 


_-m/2+l _-m/2+l 

a f _ a i 


-m/2 + 1 


(3) 


( 5 ) 


or, in general, 

N = f(Xj), i = 1 6, 

where the Xi are the six independent variables in equations (3) and (4). Equation (3) is 
"randomized" by assuming the first four variables in equation (4) to be random. Assuming 
a small crack in a relatively large component leads to assuming Y = 1.0, a deterministic 
value. A deterministic final crack size was chosen since experimental evidence indicated 
that it was relatively unimportant 1 

Probabilistic analysis, via simulation, yields the distribution of the dependent 
random variable, cycles, N. A probability density function (p.d.f.) of cycles is generated 
using the maximum penalized likelihood method. Maximum penalized likelihood generates 
the p.d.f. estimate using the method of maximum likelihood together with a penalty 
function to smooth it. 3 

2.2 RANDOM2 Computer Program 

A FORTRAN computer program for the fatigue crack growth model, called 
RANDOM2, was written using the above-described probabilistic methodology and the 
constitutive relationship expressed in equation (3). Although the four independent random 
variables could have any distribution, this initial program provided for normal or lognormal 
only. 


A complete Users Manual for RANDOM2 is contained in Appendix 1. Also, a disk 
containing a new version of RANDOM2 and a sample problem accompanies this report. 
The new version of RANDOM2, documented in the Users Manual, uses the new ISML, 
Ver. 10 subroutines and provides for parameter input from an input file. 
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3.0 FATIGUE STRENGTH REDUCTION MODEL 


3.1 Background 

Fatigue strength data are usually presented as cycles to failure for each of several 
stress amplitudes, the familiar S-N diagram. Results indicate that for lower stress 
amplitudes the cycles (or time) to failure increases. Thus, a power curve fit through the 
data yields a monotonically decreasing curve. In general, this curve is represented as 

S = [N/C'f 1/m (6) 

where the primitive variables in this equation are as follows: S is the applied constant 
amplitude alternating stress at failure or fatigue strength, N is number of cycles, C' is a 
material parameter that varies from specimen to specimen and m' is a material constant. 4 
Equation (6) can be written in terms of "cycles to reach a given fatigue strength" as 

N = C'S'" 1 ’ (7) 

Recently another fatigue strength reduction model has been proposed that takes into 
account the effect of temperature as well as other parameters that affect strength. 5 The 
general form of the constitutive relationships for this model is applied to the constituents of 
high temperature composite materials. Specifically, it is applied herein for the case of a 
single material constituent The mechanical property of interest is fatigue strength which is 
expressed in terms of primitive variables, including the general categories of temperature, 
mechanical cycles and mean stress. For these catagories, the relationship becomes 
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where S is the applied constant amplitude alternating stress at failure (fatigue strength) at 

current (or operating) temperature, T, mean stress, o, and mechanical cycle, Nm- So is 
fatigue strength at reference temperature. To (usually room temperature), reference mean 
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or melting temperature of the material, Sf is the final or tensile strength of the material, and 
Nmf is the final mechanical cycle or lifetime. Empirical parameters, n, m, and q, are 
determined from available experimental data or estimated from anticipated behavior of the 
particular product term 6 . Note that the term containing mechanical cycles is expressed in 
terms of the log of cycles rather than cycles. This formulation is attractive when Nm and 
Nmo are small compared to Nmf- The equation may be solved for Nm, or the "cycles to 
reach a given fatigue strength." The expression is 
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For values typical of a cast nickel base-superalloy subjected to typical loads and 
temperatures, equation (9) indicates increasing life for decreasing temperature, decreasing 
tensile mean stress, and decreasing applied alternating stress. It indicates decreasing life 
for increasing temperature, decreasing compressive mean stress, and increasing applied 
alternating stress. Therefore, equation (9) predicts observed trends in general. 

Probabilistic analysis, via simulation, yields the distribution of the dependent 
random variable, cycles, N. A probability density function (p.d.f.) of cycles is generated 
using the maximum penalized likelihood method for RANDOM3. For RANDOM4, a 
p.d.f. of cycles is generated using the maximum entropy method. Maximum entropy uses 
Jaynes' principle which says that "the minimally prejudiced distribution is that which 
maximizes the entropy subjected to the constraints supplied by the given information." 7 

3.2 RANDOM3 and RANDOM4 Computer Programs 

FORTRAN computer programs for the fatigue strenth reduction model called 
RANDOM3 and RANDOM4 were written using the above-described probabilistic 
methodology and the constitutive relationship expressed in equation (9). Although the 
thirteen independent random variables could have any distribution, these programs 
provided for normal or lognormal only. 

A complete Users Manual for RANDOM3 and RANDOM4 is contained in 
Appendix 2. Also, a disk containing new versions of RANDOM3 and RANDOM4 uses 
the new IMSL, Ver. 10 subroutines and provides for parameter input from an input file. 

3.3 Sensitivity Study 

The fatigue strength degradation model using the maximum entropy method of 
p.d.f. generation (RANDOM4) was selected for use in a sensitivity study. A base line 

problem utilizing a high mean stress (o = 90 ksi) and a low constant amplitude alternating 

stress at failure (S = 22.5 ksi) was established. A room temperature (T = 68°F) problem 
was executed. The input for this problem is given in Table 1 and the output, in the form of 
a p.d.f. and a c.d.f. is given in Figures 1 and 2. A high temperature (T= 1562°F) base line 
problem was also selected. Then, for a fixed base line alternating stress at failure, the mean 
stress was varied above and below the base line value. Both room and high temperatures 
were selected. Finally, for a fixed base line mean stress, the alternating stress at failure 
was varied above and below the base line value. Again, room and high temperatures were 
selected. A summary of the cases studied is given in Table 2. 

Conclusions drawn from this sensitivity study are summarized below. Increasing 
temperature for the same stress conditions reduces lifetime for all cases (see, for example, 
Figure 3). At room temperature, when mean stress is increased by 10%, lifetime decreases 
only very slightly. At high temperature, however, when mean stress is increased by 10%, 
lifetime decreases substantially. Also, at room temperature, when alternating stress is 
increased by 30%, lifetime decreases only slightly. At high temperatures, however, when 
alternating stress is increased by 30%, lifetime decreases very substantially. Considering 
the above points, lifetime is more sensitive to increasing alternating stress, rather than mean 
stress. This is probably because alternating stress was increased by 30%, whereas mean 
stress was increased by only 10%. 
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Table 1 Base iine room temperature (RT) problem input, using the fatigue strength 
reduction model with maximum entropy p.d.f. generation (RANDOM4) 

Variable Distribution Mean Std. Dev. 

Type Value % of Mean 


Tp (Melting Temp.) 

Normal 

2732.0 °F 

82.0 

3 

Sp (Ult. Tensile Str.) 

Lognormal 

130.0 ksi 

6.5 

5 

Nmf (Log of Final Cycle) 

Lognormal 

8.0 

0.8 

10 

T 0 (Ref. Temp.) 

Normal 

68.0 °F 

2.0 

3 

o 0 (Residual Comp. Stress) 

Lognormal 

-2.9 ksi 

-0.145 

5 

Nm (Log of Ref. Cycle) 

Lognormal 

7.0 

0.7 

10 

S 0 (Ref. Fatigue Str.) 

I^ognormal 

72.6 ksi 

3.6 

5 

T (Current Temp.) 

Normal 

68.0 OF 

2.0 

3 

o (Current Mean Stress) 

Lognormal 

90.0 ksi 

4.5 

5 

S (Current Fatigue Str.) 

Lognormal 

22.5 ksi 

1.125 

5 

n (Temp. Exponent) 

Normal 

0.5 

0.015 

0.3 

m (Stress Exponent) 

Normal 

0.5 

0.015 

0.3 

q (Cycle Exponent) 

Normal 

0.5 

0.015 

0.3 



log of cycles 


Fig. 1 p.d.f. of base line room temperature (RT) problem 
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log of cycles 

Fig. 2 c.d.f. of base line room temperature (RT) problem 


Table 2 Sensitivity study cases for fatigue strength reduction model 
with maximum entropy p.d.f. generation (RANDOM4) 


0 (ksi) 

S (ksi) 

T (°F) 

80 

22.5 

68 


Base Line (RT) 

90 

22.5 

68 

100 

22.5 

68 


80 

22.5 

1562 

Base Line (HT) 

90 

22.5 

1562 


100 

22.5 

1562 


90 

15.0 

68 

Base Line (RT) 

90 

22.5 

68 


90 

30.0 

68 


90 

15.0 

1562 

Base Line (HT) 

90 

22.5 

1562 


90 

30.0 

1562 
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log of cycles 


Fig. 3 c.d.f. of base line room temperature (RT) problem compared with 
c.d.f. of base line high temperature (HT) problem. 
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1.0 INTRODUCTION 


This User Manual documents the FORTRAN program RANDOM2. RANDOM2 is 
based on fracture mechanics using a probabilistic fatigue crack growth model. It predicts 
the random lifetime of an engine component to reach a given crack size (see Section 2.0, 
Theoretical Background). 

Included in this Manual are details regarding the theoretical background of 
RANDOM2, input data instructions and a sample problem illustrating the use of 
RANDOM2. Appendix A gives information on the physical quantities, their symbols, 
FORTRAN names, and both SI and U.S. Customary units. Appendix B includes 
photocopies of the actual computer printout corresponding to the sample problem. 
Appendices C and D detail the IMSL, Ver. 10 1 , subroutines and functions called by 
RANDOM2 and a S AS/GRAPH 2 program that can be used to plot both the probability 
density function (p.d.f.) and the cumulative distribution function (c.d.f.). 
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2.0 THEORETICAL BACKGROUND 


Fatigue crack growth data are usually presented as cycles, N, to reach a particular 
crack length, a. The initial crack size is a|. It is generally accepted that under constant 
amplitude alternating stress, fatigue crack growth can be related to stress intensity through a 
first order differential equation.-* 

da/dN = C(AK) m (1) 

where C is a material parameter, m is a material property (often a constant) and AK is the 
stress intensity range. Stress intensity range is given by 

AK = YAo/ia 

where Y is a constant dependent upon component and crack geometry and Aa is the 
constant amplitude alternating stress. Therefore, equation (1) can be written as 

da/dN = C(YAo/ia ) m 


or, 


da/dN = C Y m Ao m 7t m/2 a m/2 . (2) 

Equation (2) can be integrated, from the initial crack length, ai, to the final crack 
length, af, to yield N, the number of cycles. The result is 


N = 


1 

CY m Jt m/2 Ao m . 


_-m/2+l -m/2+1 

a f ' E i 


-m/2 + 1 


(3) 


Thus, equation (3) gives the "cycles to reach a given crack length." 

Metallurgical evidence indicates that casting pores play a significant role in the high- 
cycle fatigue life of cast nickel base-superalloys, especially at high temperatures. 4 The 
location and size of these fatigue crack-initiating pores vary greatly from one aerospace 
propulsion system component to another. This accounts for the large variability in fatigue 
life and leads to consideration of fatigue crack growth as a random phenomenon. 

Fatigue life directly relates to casting pore size, and pore size can be used to 
determine initial crack size, ai. Thus, utilizing principles of both probabilistic analysis and 
fatigue crack growth, a quantitative probabilistic constitutive relationship between fatigue 
life and fracture mechanics parameters can be developed. Using the "randomized equation" 
approach, the fatigue crack growth model, given by equation (3) has the following form: 

N = f(C,m,Aa,a i ,a f ,Y) 


(4) 


or, in general, 


N = f(X.), i = 1,....,6, (5) 

where the Xi arc the six independent variables in equations (3) and (4). Equation (3) is 
"randomized" by assuming the first four variables in equation (4) to be random. Assuming 
a small crack in a relatively large component leads to assuming Y = 1.0, a deterministic 
value. A deterministic final crack size was chosen since experimental evidence indicated 
that it was relatively unimportant 5 

Probabilistic analysis, via simulation, yields the distribution of the dependent 
random variable, cycles, N. A probability density function (p.d.f.) of cycles is generated 
using the maximum penalized likelihood method. Maximum penalized likelihood generates 
the p.d.f. estimate using the method of maximum likelihood together with a penalty 
function to smooth it. 5 
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3.0 INPUT DATA 


Data input for RANDOM2 is user friendly and easy to manipulate (see, for 
example, the file entitled NORMAL.INP, in Section 4.0). Die first five lines of input have 
the same format, namely 2E12.4, and the last two lines differ. The last two lines of input 
have the formats I3,2X,I3,2X,2E12.4,2X,I3 and 13, respectively. A brief line by line 
description is given along with an example for each line (Note: the ruler is to aid the user in 
formatting and is not a part of the input). A table listing the physical quantities, their units 
and symbols is given in Appendix A. 

1. Random Number Generator Seed, ISEED, and Sample Size, NTOT 
EXAMPLE: 


123456789012345678901234567890 



6.2E+02 6.2E+01 
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6. The DESPL 1 parameters are NODE, INIT, ALPHA, EPS, MAXIT and are entered 
in that order as follows: 

EXAMPLE: 

1234567890123456789012345678901234567890 
21 0 50.0E-01 10.0E-05 30 

7. The DESPL parameter, IOPT, is entered as follows: 

EXAMPLE: 

1234567890 
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4.0 SAMPLE PROBLEM FOR RANDOM2 


The objective of tins program is to predict the random lifetime, to reach a given 
crack size for an engine component . The theory is based on fracture mechanics, using a 
probabilistic fatigue crack growth model (see Section 2.0, Theoretical Background). 
RANDOM2 input parameters are given in Table Al.l. Note that the first four parameters 
are random. Their means and standard deviations are input by the user. The last two 
parameters, Af and Y, are deterministic and are fixed internally by the program. They are 
equal to the values shown in Table ALL 


Table Al.l RANDOM2 sample problem input (SI units) 


FORTRAN 

Name 

Distribution 

Type 

Mean 

Standard Deviation 
(Value) (% of Mean) 

RMM 

normal 

28.0E-01 

1.4E-01 (5%) 

AI 

lognormal 

300.0E-06 

45.0E-06 (15%) 

RCC 

lognormal 

2.20E-11 

0.22E-11 (10%) 

DELSIG 

lognormal 

6.2E+02 

6.2E+01 (10%) 

AF 

N/A 

2.0E-03 

N/A 

YY 

N/A 

1.0 

N/A 


The input is entered in the following format in a file entitled NORMAL.INP. 

1234567890123456789012345678901234567890 
1 

28.0E-01 
300.0E-06 
2 . 20E-11 
6.2E+02 
21 0 
2 


40 

1.4E-01 

45.0E-06 

0.22E-11 

6.2E+01 

50.0E-01 10.0E-05 30 


1 8 


Execution of RANDOM2 (source code entitled NR2.FOR) produces an output file 
entitled RANDM22 giving intermediate results (see Appendix B). Execution also produces 
the plotfiles OUT1 and OUT2 (see Appendix B). These files are used to plot the X and Y 
axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.d.f.), respectively, generated by RANDOM2. The plots are drawn from the plotfiles by 
the SAS/GRAPH graphing program (see Appendix C). These plots for the sample 
problem are shown in Figures A 1.1 and A 1.2. 

This same sample problem has been reported in Boyce and Chamis. 6 There, 
however, it utilized U.S. Customary units and an older version of RANDOM2 
(IMSL Version 9.2 subroutines). 



Fig. A 1 . 1 p.d.f. of log of mechanical cycles for fatigue crack 
growth model, using maximum penalized likelihood. 
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CUMULATIVE DISTRIBUTION FUNCTION 



Fig. A 1 .2 c.d.f. of log of mechanical cycles for fatigue crack 
growth model, using maximum penalized likelihood. 


| 
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6.0 APPENDIX A 


PHYSICAL QUANTITIES, SYMBOLS, AND UNITS 

The physical quantities, their symbols, and units for the fatigue crack growth 
model are given in the following table. 


Table A 1.2 Physical quantities, symbols, and units 
for fatigue crack growth model for RANDOM2 


Physical 

Theory 

FORTRAN 

Units 

Quantity 

Symbol 

Name 

SI 

U.S. 

Material Property 

m 

RMM 

m/cycle/M Pa m 

in/cycle/ksi in 

Initial Crack Size 

Ai 

RAI 

m 

in 

Material Properly 

C 

RCC 

m/cycle 

in/cycle 

Alternating Stress 

Aa 

DELSIG 

MPa 

ksi 

Final Crack Size 

Af 

AF 

m 

in 

Geometry Dependent 

Y 

YY 

(dimensionless) 


Constant 
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7.0 APPENDIX B 


SAMPLE PROBLEM: SOURCE, INPUT AND OUTPUT FILES 
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8.0 APPENDIX C 


1. RNSET 

2. RNNOR 

3. RNLNL - 

4. DESPL - 

5. GCDF - 


IMSL SUBROUTINE CALLS FROM RANDOM2 


Initializes a randan seed for use in the IMSL random number generators. 

Generates pseudorandom numbers from a standard normal distribution 
using an inverse CDF method. 

Generates pseudorandom numbers from a lognormal distribution. 

Performs nonparametric probability density function estimation by the 
penalized likelihood method. 

Evaluates a general continuous cumulative distribution function given 
ordinates of die density. 
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9.0 APPENDIX D 


SAMPLE SAS/GRAPH (VER. 5.16) PROGRAM FOR RANDOM2 


data a; 

INFILE 'OUT1.CPR' FIRSTOBS=2;input x y; 

GOPTIONS DEVICE=HP7470; 
proc gplot; 

axisl label=(h=l f=simplex 'LOG OF CYCLES’) 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A=90 F=SIMPLEX ’PROBABILITY DENSITY FUNCTION’; 
symbol i=spline v=square; 
dataB; 

INFILE ’OUT2.CPR’ FIRSTOBS=2;input x y; 
proc gplot; 

axisl label=(h=l f=simplex ’LOG OF CYCLES’) 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A =90 F=SIMPLEX ’CUMULATIVE DISTRIBUTION FUNCTION’; 
symbol i=spline v=squarc; 
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6.0 APPENDIX 2 


FATIGUE STRENGTH DEGRADATION MODEL: 
RANDOM3 AND RANDOM4 USER MANUAL 


FATIGUE STRENGTH REDUCTION MODEL: 
RANDOM3 and RANDOM4 USER MANUAL 


Prepared by : 

Lola Boyce, Ph.D., P.E. 
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APPENDIX 2 
of Annual Report 
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1.0 INTRODUCTION 


This User Manual documents the FORTRAN programs RANDOM3 and 
RANDOM4. They are based on fatigue strength reduction, using a probabilistic 
constitutive model. They predict the random lifetime of an engine component to reach a 
given fatigue strength (see Section 2.0, Theoretical Background). 

Included in this Manual are details regarding the theoretical backgrounds of 
RANDOM3 and RANDOM4, input data instructions and sample problems illustrating the 
use of RANDOM3 and RANDOM4 . Appendix A gives information on the physical 
quantities, their symbols, FORTRAN names and both SI and U.S. Customary units. 
Appendix B and C include photocopies of the actual computer printout corresponding to the 
sample problems. Appendices D and E detail the IMSL, Version 10 l , subroutines and 
functions called by RANDOM3 and RANDOM4 and SAS/GRAPH 2 programs that can be 
used to plot both the probability density functions (p.d.f.) and the cumulative distribution 
functions (c.d.f.). 
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2.0 THEORETICAL BACKGROUND 


Fatigue strength data are usually presented as cycles to failure for each of several 
stress amplitudes, the familiar S-N diagram. Results indicate that for lower stress 
amplitudes the cycles (or time) to failure increases. Thus, a power curve fit through the 
data yields a monotonically decreasing curve. In general, this curve is represented as 

S = [N/CT 1/m (6) 

where the primitive variables in this equation are as follows: S is the applied constant 
amplitude alternating stress at failure or fatigue strength, N is number of cycles, C' is a 
material parameter that varies from specimen to specimen and m' is a material constant. 3 
Equation (6) can be written in terms of "cycles to reach a given fatigue strength" as 

N = C'S' m ’ (7) 

Recently another fatigue strength reduction model has been proposed that takes into 
account the effect of temperature as well as other parameters that affect strength. 4 The 
general form of the constitutive relationships for this model is applied to the constituents of 
high temperature composite materials. Specifically, it is applied herein for the case of a 
single material constituent The mechanical property of interest is fatigue strength which is 
expressed in terms of primitive variables, including the general categories of temperature, 
mechanical cycles and mean stress. For these categories, the relationship becomes 

s r vit r V°f r [ °g n mf - '° s n m f 

So -LV T oJ L S F'°»J L 1<>SN MF" 108 N MoJ ® 

where S is the applied constant amplitude alternating stress at failure (fatigue strength) at 
current (or operating) temperature, T, mean stress, a, and mechanical cycle, Nm. So is 
fatigue strength at reference temperature, To (usually room temperature), reference mean 

stress (or residual stress), Co, and reference mechanical cycle, Nmo- Also, Tp is the final 
or melting temperature of the material, Sp is the final or tensile strength of the material, and 
Nmf is the find mechanical cycle or lifetime. Empirical parameters, n, m, and q, are 
determined from available experimental data or estimated from anticipated behavior of the 
particular product term. 5 Note that the term containing mechanical cycles is expressed in 
terms of the log of cycles rather than cycles. This formulation is attractive when Nm and 
Nmo are small compared to Nmf- The equation may be solved for Nm, or the "cycles to 
reach a given fatigue strength." The expression is 



5 1 


For values typical of a cast nickel base-superalloy subjected to typical loads and 
temperatures, equation (9) indicates increasing life for decreasing temperature, decreasing 
tensile mean stress, and decreasing applied alternating stress. It indicates decreasing life 
for increasing temperature, decreasing compressive mean stress, and increasing applied 
alternating stress. Therefore, equation (9) predicts observed trends in general. 

Probabilistic analysis, via simulation, yields the distribution of the dependent 
random variable, cycles, N. A probability density function (p.d.f.) of cycles is generated 
using the maximum penalized likelihood method for RANDOM3. For RANDOM4, a 
p.d.f. of cycles is generated using the maximum entropy method. Maximum entropy uses 
Jaynes' principle which says that "the minimally prejudiced distribution is that which 
maximizes die entropy subjected to the constraints supplied by the given information." 6 
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3.0 INPUT DATA 


Data input for RANDOM3 and RANDOM4 is user friendly and easy to manipulate 
(see, for example, the file entitled NORMAL.INP, in Section 4.0). The first twelve lines 
of input have die same format, 2E12.4 and the last two lines differ. The last two lines of 
input have the formats I3,2X,I3,2X,2E12.4,2X,I3 and 13, respectively. A brief, line by 
line description is given along with an example for each line (NOTE: the ruler is to aid the 
user in formatting and is not a part of the input). A table listing the physical quantities, 
their units and symbols is given in Appendix A. 

1. Random Number Generator Seed, ISEED, and Sample Size, NTOT 

EXAMPLE: 

123456789012345678901234567890 
1 40 


2. Ultimate Tensile Strength, SF 
EXAMPLE: 

123456789012345678901234567890 

900.0000 45.0000 


3. Log of Final Cycle, NMF 
EXAMPLE: 

123456789012345678901234567890 

8.0000 0.8000 


4. Reference Fatigue Strength, SO 
EXAMPLE: 

123456789012345678901234567890 

500.0000 25.0000 


5. Log of Reference Cycle, NMO 
EXAMPLE: 

123456789012345678901234567890 

7.0000 0.7000 
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6. Current Fatigue Strength, S 
EXAMPLE: 

123456789012345678901234567890 

250.0000 12.0000 

7. Residual Compressive Stress, SIGO 
EXAMPLE: 

123456789012345678901234567890 
20.0000 1.0000 


8. Current Mean Stress, SIG 
EXAMPLE: 

123456789012345678901234567890 

150.0000 7.5000 


9. Temperature Exponent, XXN, Stress Exponent, XXM, and Cycle Exponent, XXQ 
EXAMPLE: 

123456789012345678901234567890 
0.5000 0.0150 

10. Melting Temperature, TF 
EXAMPLE: 

123456789012345678901234567890 

1500.0000 75.0000 

11. Reference Temperature, TO 
EXAMPLE: 

123456789012345678901234567890 

20.0000 0.6000 
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12. Current Temperature, T 
EXAMPLE: 

123456789012345678901234567890 
850.0000 25.0000 


13. The DESPL 1 parameters are NODE, INIT, ALPHA, EPS, and MAXIT and are 
entered in that order as follows: 

EXAMPLE: 

1234567890123456789012345678901234567890 
21 0 20.0000 1.0E-05 30 

14. The DESPL parameter, IOPT, is entered as follows: 

EXAMPLE: 

1234567890 

2 
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4.0 SAMPLE PROBLEMS FOR RANDOM3 AND RANDOM4 


The objective of these programs is to predict the random lifetime to reach a given 
fatigue strength for an engine component. The theory is based on fatigue strength 
reduction, using a probabilistic constitutive model. The only difference between 
RANDOM3 and RANDOM4 is the method used to generate p.d.f. estimates. RANDOM3 
uses maximum penalized likelihood, while RANDOM4 uses maximum entropy (see 
Section 2.0, Theoretical Background). RANDOM 3 and RANDOM4 input parameters are 
given in Table A2.1. 


TABLE A2. 1 RANDOM3 and RANDOM4 input (SI units) 


FORTRAN 

Name 

Distribution 

Type 

Mean 

Standard Deviation 
(V alue) (% of Mean) 

SF 

normal 

900.0 

45.0 

(3%) 

NMF 

lognormal 

8.0 

0.8 

(10%) 

SO 

lognormal 

500.0 

25.0 

(5%) 

NMO 

lognormal 

7.0 

0.7 

(10%) 

S 

lognormal 

250.0 

12.5 

(5%) 


lognormal 

-20.0 

-1.0 

(1%) 


lognormal 

150.0 

7.5 

(5%) 


normal 

0.5 

0.015 

(0.3%) 

XXM 

normal 

0.5 

0.015 

(0.3%) 

XXQ 

normal 

0.5 

0.015 

(0.3%) 

TF 

normal 

1500.0 

45.0 

(3%) 

TO 

normal 

20.0 

0.6 

(3%) 

T 

normal 

850.0 

25.5 

. (3%) 
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The input is entered in the following format in a file 
entitled NORMAL.INP. 

1234567890123456789012345678901234567890 


1 

■HKIi] 

900.0000 

45.0000 

8.0000 

0.8000 

500.0000 

25.0000 

7.0000 

0.7000 

250.0000 

12.5000 

20.0000 

1.0000 

150.0000 

7.5000 

0.5000 

0.0150 

1500.0000 

75.0000 

20.0000 

0.6000 

850.0000 

25.5000 

0 

20.00 


2 


Execution of RANDOM3 and RANDOM4 (source code entitled NR3.FOR and 
NR4.FOR, respectively) produces files entitled RANDM33 and RANDM44. These give 
intermediate results (see Appendices B and C). Execution also produces plotfiles entitled 
PLOT1 and PLOT2 (see Appendices B and C). These files are used to plot the X and Y 
axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.cLf.), respectively, generated by RANDOM3 and RANDOM4. The plots are drawn 
from tiie plotfiles by the SAS/GRAPH graphing program (see Appendix D). These plots 
for the sample problem are shown Figures 1, 2, 3, and 4. This same sample problem has 
been reported in Boyce and Chamis. 7 There, however, it utilized U.S. Customary units 
and older versions of RANDOM3 and RANDOM4 (using IMSL Version 9.2 subroutines) 
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I 

I 

I 


' ' I I 

3 8 7 8 9 10 

LOG OF CYCLES 

Fig. A2.2 c.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum penalized likelihood method of p.d.f. generation. 
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Fig. A2.3 p.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum entropy method of p.d.f. generation. 



Fig. A2.4 c.d.f. of log of mechanical cycles for fatigue strength reduction model, 
using maximum entropy method of p.d.f. generation. 
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6.0 APPENDIX A 


PHYSICAL QUANTITIES, SYMBOLS, AND UNITS 


The physical quantities, their symbols and units for the fatigue crack growth 
model are given in the following table. 


Table A2.2 Physical quantities, symbols, and units for fatigue 
crack growth model for RANDOM3 and RANDOM4. 


Physical 

Quantity 

Theory 

Symbol 

FORTRAN 

Name 

Units 

SI U.S 

Ultimate Tensile 





Strength 

SF 

SF 

MPa 

ksi 

Final Cycle 





(lifetime) 

n mf 

NMF 

dimensionless 

Reference Fatigue 





Strength 

SO 

SO 

MPa 

ksi 

Reference Cycles 

n mo 

NMO 

dimensionless 

Current Fatigue 





Strengths 

S 

S 

MPa 

ksi 

Residual Compressive 





Stress 

Oo 

SIGO 

MPa 

ksi 

Current Mean Stress 

a 

SIG 

MPa 

ksi 


n 

XXN 

dimensionless 

Empirical Material 





Parameters 

m 

XXM 

dimensionless 


q 

XXQ 

dimensionless 

Melting Temperature 

TF 

TF 

°C 

°F 

Reference 





Temperature 

TO 

TO 

°C 

°F 

Current Temperature 

T 

T 

°C 

°F 
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7.0 APPENDIX B 


RANDOM3 SAMPLE PROBLEM: SOURCE, INPUT AND OUTPUT FILES 


! 
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c w 


DUO 
DUO 
(j'jij 
DUO 
DO O 

ouu 
DO fj 
DUO 


f ir in 
rK'fj 
r mcmc j 
t ir k j 
cm fK j 
cicm cm 
r Mr Mf i 
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i ,*r Mr m 
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CJOO 

U-’U 


Dm 

mm 

M D 

mm 

mm 
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U> 
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mm 

mm 


mm 

mm 


mm 

mm 

mm 

mm 

mm 
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UJ 
UJ 
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UJ I iJ 
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Ul UJ 
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mm 
mm 
m m 
mm 
mm 
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mmmm 


mmmm 

mmmm 


mm 

mm 

mm 

mm 

mm 

mm 


mm 

mm 


mm 

mm 


2 

2 

2 

2222222 


CJ 

u 

o 


u 

CJ 

u 


u DU 0*0 


UJ 

LU 


LU 

LU LU 

UJ LU 

UJ UJ UJ 

LU Ul LU LU LU LU LU 
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ct et ct o' lycy 
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CY CY U- uY lY CY 
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cc o cycy 
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22222222222222 

22222222222222 

22 

22 

22 

22 

22 

22 

22222222222222 


cycycy cy lY a: 

CYCYlYCY occy 

CY LY CYCY CYCY 

cycy cyo cycy 

a: a: zcczcc 
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lYCY CYCY 

cycy ecu: 

cycy cycy cycy cycy cycy cycy cycy 

cy ec cc cy ec cc ec cc cc cy oc cy ec cc 


OOOOOOOOOO 
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oo oo 

oo oo 
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oo oo 

OO OO 
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RANDOM4 SAMPLE PROBLEM: SOURCE, INPUT AND OUTPUTFTLES 
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9.0 APPENDIX D 

IMSL SUBROUTINE CALLS FROM RANDOM3 AND RANDOM4 


RANDOM3 

1. RNSET - 

2. RNNOR - 

3. RNLNL - 

4. DESPL - 

5. GCDF - 

RANDOM4 

1. RNSET 

2. RNNOR 

3. RNLNL 


Initializes a random seed for use in the IMSL 
random number generators. 

Generates pseudorandom numbers from a standard 
normal distribution using an inverse CDF 
method. 

Generates pseudorandom numbers from a 
lognormal distribution. 

Performs nonparametric probability density 
function estimation by the penalized 
likelihood method. 

Evaluates a general continuous cumulative 
distribution function given the ordinates 
of the density. 


Initializes a random seed for use in the IMSL 
random number generators. 

Generates pseudorandom numbers from a standard 
normal distribution using an inverse CDF 
method. 

Generates pseudorandom numbers from a 
lognormal distribution. 
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10.0 APPENDIX E 


SAMPLE SAS/GRAPH PROGRAM FOR RANDOM3 AND RANDOM4 


data a; 

INFILE 'PLOT1.CPR' FIRSTOBS=2;input x y; 

GOPTIONS DE VICE=HP7 47 0; 
proc gplot; 

axisl label=(h=l f=simplex ’LOG OF CYCLES’) 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A =90 F=SIMPLEX ’PROBABILITY DENSITY FUNCTION'; 
symbol i=spline v=square; 
dataB; 

INFILE ’PLOT2.CPR’ FIRSTOBS=2;input x y; 
proc gplot; 

axisl label=(h=l f=simplex 'LOG OF CYCLES') 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A =90 F=SIMPLEX ’CUMULATIVE DISTRIBUTION FUNCTION'; 
symbol i=spline v=square; 
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